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The corrosion of steel in concrete is mainly due to the chemical reaction between the chloride ions and iron 
ions. Indeed, this is a serious threaten for reinforced concrete structure, especially for the reinforced 
concrete structure in the sea. So it is urgent and important to protect concrete against chloride ions 
corrosion. In this work, we report multilayer concrete can cloak chloride ions. We formulated five kinds of 
concrete A, B, C, D and E, which are made of different proportion of cement, sand and glue, and fabricated 
six-layer (ABACAD) cylinder diffusion cloak and background media E. The simulation results show that the 
six-layer mass diffusion cloak can protect concrete against chloride ions penetration, while the experiment 
results show that the concentration gradients are parallel and equal outside the outer circle in the diffusion 
flux lines, the iso-concentration lines are parallel outside the outer circle, and the concentration gradients in 
the inner circle are smaller than those outside the outer circle. 

I n 2006, Pendry designed perfect invisibility cloak 1 . Following this method, varieties of electromagnetic cloak 2-17 , 
I acoustic cloak 18 23 , mater wave cloak 24 , magnetic cloak 25 27 , dc electric cloak 28 ' 29 , heat cloak 30-38 and mass dif- 
I fusion cloak 39 have been designed. The first experiment 4 demonstrated invisible cloak having adopted meta- 
material with simplified material parameters, the simplified cloak inherits some properties of the ideal cloak, but 
finite scatterings exist. Perfect invisible cloak requires inhomogeneous and anisotropic media, however, it is 
difficult to construct by using natural materials. Recently, a cylinder heat cloak structure that does not require 
metamaterials to realize the anisotropy and inhomogeneity of the material parameters was proposed 32 ; it was 
realized by a concentric layered structure consisting of alternating homogeneous isotropic materials, which can be 
treated as an effective medium with the radius dependent anisotropy. The time-dependent thermodynamic 
problem has also been treated theoretically 33 and experimentally 37 by Guenneau and co-author. Subsequent 
Guenneau 39 adapted the concept of transformational thermodynamics, the n-dimensional, time- dependent, aniso- 
tropic heterogeneous Ficks equation is derived. Based on the above research, we fabricated six-layer cylinder mass 
diffusion cloak according to the diffusion theory. The experiment demonstrated that the six-layer structures can 
also cloak diffusion flux. Our work proves the feasibility to extend the electromagnetic transformation optic devices 
into the mass diffusion area, which may pave a way to control the diffusion flux for novel device applications. 

Results 

Mass diffusion principle. The mass diffusion equation can be written as: 

dC 



dt 



= V-(DVQ 



(1) 



where t is time, C is mass diffusion concentration (mol • m 3 or kg • m 3 , et. al), D is diffusion coefficient (m 2 • s x ), 

V=-T+-7*+ - k is the operator in the Cartesian coordinate system. Generally speaking, the mass diffusion 

x y z 

concentration is the function of space and time coordinate (x,y, z, t), i.e., C = C (x,y, z, t). According to the mass 
diffusion theory, we can obtain the following simplified radial and tangential mass diffusion coefficient similar 
thermal diffusion coefficient 33,37 
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Figure 1 | Principle of mass diffusion, (a) The relative diffusion coefficient of the hollow cylinder (D c = 1.0) is equal to the relative diffusion coefficient of 
the background media (D b = 1.0). (b) The hollow cylinder with large relative diffusion coefficient (D c = 4.0) attracts diffusion flux lines, (c) The hollow 
cylinder with giant relative diffusion coefficient (D c = 10000) completely attracts diffusion flux lines, (d) The hollow cylinder with small relative 
diffusion coefficient (D c = 0.25) repels diffusion flux lines, (e) The hollow cylinder with very small relative diffusion coefficient (D c = 0.000 1 ) completely 
repels diffusion flux lines, (f) Mass diffusion cloak, the diffusion coefficient of the hollow cylinder fulfills Eq. (2). 



The cylindrical anisotropy can be obtained by layered structure of 
homogeneous isotropic materials 33 . Assumed that two kinds of 
materials have a homogeneous and an isotropic diffusion coeffici- 
ents D A and D B , respectively, the thicknesses are d A and d B , respec- 
tively, we obtain the following effective parameters D r and D e 
d A + d B _ d A d B 

-Dr = D- A iv (3) 

_ D A d A + D B d B 
d A + d B 

where D r is the radial diffusion coefficient, and Dq is the tangential 
diffusion coefficient. 



According to the diffusion theory, we simulate mass diffusion. 
Fig. la shows that the hollow cylinder and the background medium 
have identical relative diffusion coefficient (D c = D b = 1.0), the 
diffusion coefficient of the hollow cylinder may be an effective dif- 
fusion coefficient of isotropy and homogeneous composite material, 
the diffusion flux lines are parallel both in the background media and 
hollow cylinder, this phenomenon is called "mass diffusion trans- 
parency". Fig. lb shows that the relative diffusion coefficient of the 
hollow cylinder (D c = 4.0) is larger than that of the background 
media (D b = 1.0), the diffusion flux lines are attracted by the hollow 
cylinder, and the diffusion flux lines are distorted both inside and 
outside the hollow cylinder. Fig. lc shows that the relative diffusion 




Figure 2 | Simulate steady state six-layer diffusion cloak, (a) The surface 
and the iso-concentration lines. The left is large concentration, the right is 
small concentration. The vertical lines are the iso-concentration lines, 
(b) The horizontal lines are diffusion flux lines, and the diffusion flux is 
opposite to the direction of the concentration gradient. 




Figure 3 | Simulate steady state ten-layer diffusion cloak, (a) The surface 
and the iso-concentration lines. The left is large concentration, the right is 
small concentration. The vertical lines are the iso-concentration lines, 
(b) The horizontal lines are diffusion flux lines; the performance of ten- 
layer diffusion cloak is better than that of six-layer diffusion cloak. 
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Figure 4 | Simulate time-dependent six-layer diffusion cloak, (a) The 

surface, iso-concentration lines and diffusion flux arrows. The vertical 
lines are the iso-concentration lines which are parallel outside the outer 
circle. The horizontal arrows represent diffusion fluxes which are parallel 
outside the outer circle. The left is large concentration and the right is small 
concentration. The diffusion time is 100 minutes, (b) The diffusion time is 
200 minutes. 

coefficient of the hollow cylinder (D c = 10000) is much larger than 
that of the background media (D b = 1.0), the diffusion flux lines are 
completely attracted by the hollow cylinder, no diffusion flux lines 
penetrate into the inner circle, this phenomenon is called "mass 
diffusion shielding", just like magnetostatic (electrostatic) shielding 
is generated by large permeability (permittivity) material. Fig. Id 
shows the relative diffusion coefficient of the hollow cylinder (D c 
= 0.25) is smaller than that of the background media (D b = 1.0), 
the diffusion flux lines are repelled by the hollow cylinder, and the 
diffusion flux lines are also distorted both inside and outside the 
hollow cylinder. Fig. le shows the relative diffusion coefficient of 
the hollow cylinder (D c = 0.0001) is much smaller than that of the 
background media (D b = 1.0), the diffusion flux lines are completely 
repelled by the hollow cylinder, no diffusion flux lines penetrate into 
the inner circle, just like magnetic field lines are completely repelled 
by an "superconductor". Fig. If shows the mass diffusion cloak, the 
diffusion coefficient of the hollow cylinder fulfills Eq. (2), the dif- 
fusion flux lines are parallel outside the outer cylinder, no diffusion 
flux lines penetrate into the inner circle. Theoretically and experi- 
mentally, the conditions of the mass diffusion cloak are that the 




Figure 5 | Simulate time-dependent non-cloak, (a) The surface, iso- 
concentration lines and diffusion flux arrows. The vertical lines are the iso- 
concentration lines which are not parallel outside the outer circle. The 
horizontal arrows are diffusion fluxes which are not parallel outside the 
outer circle. The diffusion time is 100 minutes, (b) The diffusion time is 
200 minutes. 

concentration gradients are parallel and equal outside the outer cir- 
cle, the iso-concentration lines are parallel outside the outer cylinder, 
the concentrations and concentration gradients in the inner circle are 
smaller than those outside the outer circle. 

Simulate steady mass diffusion cloak. According to Eq. (2) and Eq. 

(3), we design six-layer (ABAC AD) diffusion cloak and simulate 
concentration distributions. The inner and outer radii of the six- 
layer (ABACAD) are (1, 1.403), (1.403, 2), (2, 2.356) (2.356, 3), (3, 
3.318), (3.318, 4), respectively, and the diffusion coefficients of the 
four kinds of materials are D A = 4.000 X 10" 8 cm 2 /s, D B = 0.278 X 
10" 8 cm 2 /s, D c = 0.547 X 10" 8 cm 2 /s, D D = 0.741 X 10" 8 cm 2 /s, 
respectively, the diffusion coefficient of the background media are D E 
= 1.000 X 10~ 8 cm 2 /s. Fig. 2 shows the surface, iso-concentration 
lines and diffusion flux lines of the six-layer diffusion cloak; the 
vertical lines are iso-concentration lines (Fig. 2a), the red parts 
represent large concentration region, the concentration is set to be 
0.005 mol/m 3 on the left plane, and the blue parts represent small 
concentration region, the concentration is set to be 0 on the right 
plane, the horizontal lines are diffusion flux lines (Fig. 2b). The 
incident diffusion flux lines are guided around the inner circle, 
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Figure 6 | Geometry of the six-layer diffusion cloak. The diffusion cloak is in the center of the sink, the left is 0.5% NaCl solution, the right is pure water. 
The red, blue, green, black and yellow regions represent concrete A, B, C, D and E, respectively. 
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Figure 7 | Test points of the time-dependent concentrations in the iso- 
concentration lines of the six-layer cloak, (a) The parameters and the 
diffusion time are the same as that in Fig. 4a, the small blue cylinders are 
located at the iso-concentration lines, and the measurement results of the 
concentrations are shown in Fig. 8a. (b) The parameters and the diffusion 
time are the same as that in Fig. 4b, the measurement results of the 
concentrations are shown in Fig. 8b. 

outside the cloaked region, the incident diffusion flux lines return to 
their original direction as if nothing has happened. No diffusion flux 
lines penetrate into the inner circle. The diffusion flux is opposite to 
the direction of the concentration gradient. 

Fig. 3 shows the ten-layer (ABACADAEAF) diffusion cloak, 
Fig. 3a shows the surface and iso-concentration lines, Fig. 3b shows 
diffusion flux lines, the inner and outer radii of the ten-layer 
(ABACADAEAF) are (1, 1.317), (1.317, 2), (2, 2.271), (2.271, 3), 
(3, 3.236), (3.236, 4), (4, 4.208), (4.208, 5), (5, 5.186), (5.186, 6), 
respectively, and the diffusion coefficients of the six kinds of materials 
are D A = 4.000 X 1(T 8 cm 2 /s, D B = 0.253 X 1(T 8 cm 2 /s, D c = 0.478 
X 10" 8 cm 2 /s, D D = 0.650 X 10" 8 cm 2 /s, D E = 0.766 X 10" 8 cm 2 /s, 
D F = 0.855 X 10~ 8 cm 2 /s, respectively, the diffusion coefficients of 
the background media are D G = 1.000 X 10" 8 cm 2 /s. The diffusion 



Figure 9 | Test points of time-dependent concentrations in the vertical 
straight lines of non-cloak, (a) The parameters and the diffusion time are 
the same as that in Fig. 5a, the small blue cylinders are located at the vertical 
straight lines, and the measurement results of the concentrations are 
shown in Fig. 10a. (b) The parameters and the diffusion time are the same 
as that in Fig. 5b, the measurement results of the concentrations are shown 
in Fig. 10b. 

coefficient of the background media are D E = 1.000 X 10~ 8 cm 2 /s. 
Simulation results show that the more the layers, the better perform- 
ance the mass diffusion cloak, and the small diffusion coefficient 
materials exclude diffusion flux lines and the large diffusion coef- 
ficient materials attract diffusion flux lines, just like small permittivity 
(permeability) material excludes electric (magnetic) field line and 
large permittivity (permeability) material attracts electric (magnetic) 
field line. It should be pointed out that the mass diffusion cloak 
mainly acts as a flux-guiding shell to suppress the influence of dif- 
fusion reflection. 

Simulate time-dependent multilayer diffusion cloak. Utilizing the 
transient diffusion equation of COMSOL Multiphysic software, we 
simulate the diffusion cloak in different time. The parameters are the 
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Figure 8 | Measurement results of the time-dependent concentrations in the iso-concentration lines of six-layer cloak, (a) The concentrations in three 
parallel lines. The parameters and the diffusion time are the same as that in Fig. 7a. (b) The concentrations in four parallel lines. The parameters 
and the diffusion time are the same as that in Fig. 7b. 
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Figure 10 | Measurement results of the time-dependent concentrations in the vertical straight lines of non-cloak, (a) The concentrations in three 
parallel lines. The parameters and the diffusion time are the same as that in Fig. 9a. (b) The concentrations in four parallel lines. The parameters and the 
diffusion time are the same as that in Fig. 10b. 



same as that in Fig. 2. The diffusion time is 100 minutes (Fig. 4a) and 
200 minutes (Fig. 4b). Outside the outer circle, the vertical iso- 
concentration lines are parallel and the horizontal diffusion flux 
arrows are also parallel. For comparison, we also simulate the non- 
cloak in different diffusion time. The diffusion time is also 100 
minutes (Fig. 5a) and 200 minutes (Fig. 5b). The inner and outer 
radii of the hollow cylinder are a = 1, and b = 4, respectively. The 
diffusion coefficients of the background media and hollow cylinder 
are D b = 1.000 X 10" 8 cm 2 /s and Dc = 2.500 X 10" 8 cm 2 /s, 
respectively. Outside the outer circle, the iso- concentration lines 
are not parallel and the horizontal diffusion flux arrows are also 
not parallel. So this hollow cylinder is a non-cloak. 

Experiment realization. The geometry of the six-layer diffusion 
cloak is shown in Fig. 6, the geometric parameters and material 
parameters in the experiment completely accord with those in 
simulation. The red, blue, green, black and yellow regions 
represent concrete A, B, C, D and E, respectively. The diffusion 
cloak is in the center of the sink, the left is 0.5% NaCl solution, 
and the right is pure water. The concretes are made of cement, 
sand and glue. After many times try, we obtain the concretes with 
various diffusion coefficients we expect. The proportion of cement, 
sand and glue are about A (1 : 3 : 0.02), B (1 : 3 : 0.30), C (1 : 3 : 0.20), D 
(1 : 3 : 0.10), E (1 : 3 : 0.08), respectively. The diffusion coefficients are 
D A = 4.00 X 10" 8 cmVs, D B = 0.28 X 10" 8 cm 2 /s, D c = 0.55 X 
10" 8 cm 2 /s, D D = 0.74 X 10" 8 cm 2 /s, D E = 1.00 X 10" 8 cm 2 /s, 
respectively. According to the measurement results, we find that 




Figure 1 1 | Test points of the time-dependant concentration gradients in 
horizontal lines. The diffusion time is 180 minutes; the horizontal lines 
represent diffusion flux lines. The measurement results are shown in 
Fig. 12. 



the more the glue, the smaller the diffusion coefficient, i.e., the glue 
can decrease diffusion coefficient. It is also proved that other 
material, such as fly ash, can also decrease diffusion coefficient. 
According to the simplified parameter Eq. (2) and effective 
parameter Eq. (3), we fabricate the six-layer (ABACAD) cylinder 
diffusion cloak. The inner and outer radii of the six-layer cylinder 
cloak are (1 cm, 1.4 cm), (1.4 cm, 2.0 cm), (2,0 cm, 2.4 cm), 
(2.4 cm, 3.0 cm), (3.0 cm, 3.3 cm), (3.3 cm, 4.0 cm), respectively, 
the height of the cylinder cloak is h = 6 cm (in order to measure 
the mass diffusion in different time, we fabricated many identical 
diffusion cloaks which are made of the same materials). After 100 
minutes diffusion, pasting simulation results (Fig. 4a and Fig. 4b) on 
the surface of samples, we dug up a lot of cylinders (the diameter is 
4 mm) in the vertical straight lines (Fig. 7a), and measured the 
chloride ions concentration by NJCL Chloride Ions Concentra- 
tion Fast Measurement Equipment. The concentrations of 18 test 
points are shown in Fig. 8a, the concentrations in a straight line are 
nearly equal and the concentrations in three vertical straight lines 
(Fig. 8a) are about 0.0330%, 0.0280% and 0.0220%, respectively. So 
the three vertical straight lines are the iso -concentration lines. After 
200 minutes diffusion, we measured the concentrations of 24 test 
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Figure 12 | Measurement results of the time- dependant concentration 
gradients in horizontal lines. The parameters and diffusion time are the 
same as that in Fig. 11, the diffusion flux are parallel and equal. 
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Figure 13 | Measurement results of the time-dependent concentration gradients, (a) The parameters are the same as that in Fig. 4 (cloak). The red circles 
represent the concentration gradient between ( — 5.8 cm, 0) and ( — 4.2 cm, 0), and the black squares represent the concentration gradient between 
( — 0.8 cm, 0) and (0.8 cm, 0). (b) The parameters are the same as that in Fig. 5 (non-cloak). The red circles represent the concentration gradients between 
( — 5.6 cm, 0) and (—4.2 cm, 0), and the black squares represent the concentration gradients between ( — 0.8 cm, 0) and (0.8 cm, 0). 



points (Fig. 7b), the concentrations in a straight line are nearly equal 
and the concentrations in four vertical straight lines (Fig. 8b) are 
about 0.0360%, 0.0330%, 0.0180% and 0.0130%, respectively. So 
the four vertical straight lines are the iso- concentration lines. 

For comparison, we measured the chloride ions concentration in 
three parallel straight lines of non-cloak (Fig. 9a). The concentrations 
of 18 test points are shown in Fig. 10a (the measurement time is 100 
minutes), and the concentrations in a straight line are not equal. The 
concentrations of 24 test points (Fig. 9b) are shown in Fig. 10b (the 
measurement time is 200 minutes), and the concentrations in a 
straight line are not equal either, i.e., outside the outer circle, the 
iso-concentration lines are not parallel, i.e. the vertical parallel lines 
are not iso-concentration lines and the iso-concentration lines are 
not parallel. So this hollow cylinder is a non-cloak. 

By the same method, we measured the concentration gradients in 
the six horizontal lines (Fig. 11). The concentrations in the horizontal 
lines are shown in Fig. 12. The diffusion time is 180 minutes. 
Measurement results show that the concentration gradients are 
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Figure 14 | Measurement results of the time-dependent concentration in 
the center. The parameters of the red circles are the same as that in Fig. 13a 
(cloak), the parameters of the black five pointed star are the same as that in 
Fig. 13b (non-cloak). The concentrations in the center of the cloak are 
always smaller than those of the non-cloak even after 550 minutes. 



nearly equal and about 0.0027%/cm, so the horizontal lines are the 
diffusion flux lines. 

Finally, we measured time- dependent concentration gradients in 
the inner circle and outside the outer circle. We fabricated ten ident- 
ical diffusion cloak and ten identical hollow cylinders (non-cloak). In 
the inner circle, the test points are located at ( — 0.8 cm, 0) and 
(0.8 cm, 0), respectively, outside the outer circle, the test points are 
located at ( — 5.8 cm, 0) and (—4.2 cm, 0), respectively. The concen- 
tration gradients in the inner circle are always smaller than those 
outside the outer circle in the cloak even after 550 minutes (Fig. 13a). 
The concentration gradients in the inner circle are nearly equal to 
those outside the outer circle in the non-cloak (Fig. 13b). We also 
measured the time- dependent concentrations in the center of the 
cloak and non-cloak (Fig. 14). The concentrations in the center of 
the cloak are always smaller than those of the non-cloak even after 
550 minutes. The six-layer cylinder we fabricated fulfills the condi- 
tions of cloak. 

Guenneau 39 theoretically studied two and three dimensional mass 
diffusion cloak in fluids, and numerically simulated multilayer nano- 
size reduced parameter cylinder and sphere mass diffusion cloak. In 
their work, they consider the convection-diffusion equation that is a 
parabolic partial differential equation combining the diffusion equa- 
tion and the advection equation. The issue of convection is very 
importance for fluids diffusion. However, it is not necessary to con- 
sider the convection in the diffusion of chloride ions in solid con- 
crete. Other diffusion of chemical element in containers such as 
cobalt, uranium can also be cloaked by our method. The prospects 
of mass diffusion cloak are bright in the construction field. If the 
materials come into being, it would prolong the life of construction 
materials. Meanwhile it also means a lot in chemical element dif- 
fusion research. Our works provide the theory and experiment guid- 
ance to protect concrete against chloride ions corrosion. 

Discussion 

For an ideal mass diffusion cloak, the diffusion matters do not pen- 
etrate into the inner circle, thus a perfect diffusion protection is 
offered. The multilayer mass diffusion cloak is a simplified approx- 
imate cloak. The chloride ions can penetrate into the cloak region in 
long time diffusion. In a steady state, the cloak region eventually 
achieves a constant concentration. Therefore, the mass diffusion 
cloak can only provide a temporal diffusion protection rather than 
a permanent protection. Another case is that the diffusion matter 
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chloride ions pollution would change the effective parameters of the 
cloaking materials; finally the chloride ions go into the inner cylinder 
after long time diffusion. In our experiment, we chose 0.5% NaCl 
solution. The chloride ions concentration in point P (— 4.2 cm, 0) of 
Fig. 6 is less than 0.0400%, so we can ignore the influence of chloride 
ions pollution. 

Chloride ions pollution will corrode rebars. Controlling chloride 
ions diffusion in concrete has a great significance in the field of 
building material. Concrete is composed of cement, aggregate and 
paste interfacial zone 40,41 , the larger the aggregate, the smaller the 
diffusion coefficient. Increasing or decreasing diffusion coefficient 
is a traditional method for protecting concrete 42 " 44 , Actually, this 
method is a "mass diffusion shielding", just like electrostatic and 
magnetostatic shielding. We adapted the transformation optics 
method and designed six-layer mass diffusion cloak, the experiment 
demonstrated that the mass diffusion cloak can shield chloride ions 
and guide the chloride ions around the cloaked region. By mixing 
diffusion material, one can obtain various materials with a wide 
range diffusion coefficient, which is hardly available in the design 
of an electromagnet device. 

Methods 

Numerical simulation. We numerically simulate the steady and transient mass 
diffusion cloak by COMSOL Multiphysic software. The parameters fulfill Eq. (2) and 
Eq. (3). 

Fabricate samples. The samples are fabricated by six- layer concretes (ABACAD) and 
background media E, the concretes are made of cement, sand and glue mixture. The 
proportion of cement, sand and glue are A (1 : 3 : 0.02), B (1 : 3 : 0.30), C (1 : 3 : 0.20), D 
(1 : 3 : 0.10), E (1 : 3 : 0.08), respectively. 

Experiment measure. According to the simulation results of time-dependent mass 
diffusion cloak, we dug up a lot of small cylinders in the iso-concentration lines and 
concentration gradient lines, and measured the concentration of chloride ions by 
NJCL Chloride Ions Concentration Fast Measurement Equipment. 
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